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Abstract 

We explain why quantum adiabatic evolution and simulated annealing perform simi- 
larly in certain examples of searching for the minimum of a cost function of n bits. In 
these examples each bit is treated symmetrically so the cost function depends only on 
the Hamming weight of the n bits. We also give two examples, closely related to these, 
■ where the similarity breaks down in that the quantum adiabatic algorithm succeeds in 

polynomial time whereas simulated annealing requires exponential time. 

(N 
O . 

"th ■ 1 Introduction 

Quantum adiabatic evolution algorithms QXJ] are designed to minimize a (classical) cost func- 
tion whose domain is the 2 n values taken by n bits. To test the algorithm it is natural to 
look at problems where classical local search algorithms, such as simulated annealing, have 
difficulty. It is easy to construct examples where both classical local search and quantum 
adiabatic evolution require time exponential in n. An example is the Grover problem, where 
the function takes the value at a single input and is 1 on the 2 n — 1 other inputs. A variant 
of this Grover example can be found in [[|. 

More interesting examples arise when the cost function is local in the sense that it can be 
written as a sum of terms each of which involves only a few bits. A 3sat example where both 
the classical local search and "standard" quantum adiabatic algorithms are exponentially 
slow has recently been given in ||. In Section |] below we analyze an example that is based 
on H| and give a quantum mechanical argument for why quantum adiabatic algorithms have 
trouble with these kinds of examples. 

In Section [5] we give a nonlocal example where simulated annealing requires exponential 
time but quantum adiabatic evolution succeeds in polynomial time. In Section ^| we give a 
detailed analysis of a local (2s at) example that first appeared in [|IJ. Here again a quan- 
tum adiabatic algorithm succeeds even though simulated annealing fails in polynomial time. 
These examples show that there is no general theorem that quantum adiabatic algorithms 
must fail if simulated annealing fails. 

*f arhiOmit . edu, goldston@mit.edu 
^sgutmOneu. edu 
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2 Quantum adiabatic algorithms 

Quantum adiabatic evolution algorithms offer a general approach to solving combinatorial 
search problems on a quantum computer. A local combinatorial search problem can be cast 
as the classical problem of finding the minimum of cost function h(zi, . . . ,z n ) where each 
Zi — 0,1 and 

h = J2hc (1) 

c 

and each he is a nonnegative integer- valued function that depends only on a few z\. A 
quantum adiabatic algorithm converts this into the problem of producing the ground state 
of a quantum Hamiltonian Hp given by 

H P \zi, ... ,z n ) = h(z u . . . ,Zn)\zi, ... ,z n ) (2) 

where {\zi, . . . , z n ^} is a basis for the 2 n -dimensional Hilbert space of the quantum computer. 
Specifying the algorithm consists of choosing a smoothly varying Hamiltonian H (s) for < 
s < 1 such that H(l) = Hp and where -f^(O) has a ground state |^o) that is known and 
easily constructed. A run time T must also be specified. The Hamiltonian that governs the 
evolution is given by 

Hit) = H(t/T) (3) 
and the state of the quantum computer evolves according to the Schrodinger equation 

if t \m) = H{t)\m) (4) 

with the state at time t = given by 

\m) = w • (5) 

Under conditions that generally apply in the cases of interest here, the quantum adiabatic 
theorem guarantees that for T large enough the state of the quantum computer, f° r 
< t < T, will be close to the (instantaneous) ground state of H(t), and in particular |-?/>(T)^ 
will be close to the ground state of H P , encoding the solution to the problem at hand. When 
applied to a particular combinatorial search problem, the quantum algorithm is considered 
to be successful if the required running time grows polynomially in the number of bits n. 

The required running time can be related to the spectrum of H(s), in particular the 
difference between the two lowest eigenvalues E (s) and Ei(s). The required running time T 
must obey 

T » (6) 
gap 2 

where 

gap = ndn^OO - E (s)) (7) 

and £ is less than the largest eigenvalue of H P —H(0), always polynomial in n in the examples 
here. 
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3 Simulated annealing 

Simulated annealing is a classical local search strategy that can be used to find the global 
minimum of a function h(zi, . . . , z n ) of the form ([[]). For any temperature r, < r < oo, 
define the Boltzmann distribution by 

P r (zi,... ,z n ) = c(t) exp[-h(z u . . . ,z n )/r] (8) 

where c(r) is the normalizing constant such that 

£ P T (zi,... ,z n ) = l . (9) 

21,... ,Z„ 

For r = oo, all 2 n strings zi, . . . , z n are equally likely. For r = 0, the Boltzmann distribution 
concentrates on the global minimum (or minima) of n. 

The idea is to construct a Markov chain that starts in the r = oo Boltzmann distribution 
and gradually moves through the Boltzmann distributions for decreasing r down to r near 0. 
If the process succeeds, the final distribution will be close to the zero temperature Boltz- 
mann distribution, which means that the global minimum of h has been found with high 
probability. More specifically, choose a sequence of temperatures oo > t\ > r 2 > ■ ■ ■ tm > 
at which bit values may change. The r, can be chosen deterministically or randomly. Also 
choose a sequence of bits 21,22, • • • Bit if. is the bit that might be changed (according 
to an acceptance rule to be given below in (PJ)) at step k. The sequence 11,1%, ... ,1m can 
be deterministic or random. 

The initial string is picked uniformly, i.e., all 2 n strings zi, . . . ,z n have probability 
2~ n . At step k, if the current string is z±, . . . , Zi k , . . . , z n we flip bit i& and move to string 
zi, . . . , (1 — Zi k ), . . . ,z n with probability 

minjl, exp(— A/i/rfc) } 
where Ah = h(z 1 , . . . , (1 - z ik ), . . . , z n ) - h(z u . . . , z ik , . . . , z n ) . (10) 

(This is the Metropolis rule.) Otherwise the string remains unchanged. Note that bit is 
always flipped if the change decreases h, and is flipped with only small probability if the 
change would increase h by a large amount relative to the current temperature. 

The acceptance rule ( |TUD guarantees that if many transitions are made at a single tem- 
perature r, the distribution of the string will converge to the Boltzmann distribution at 
temperature r. If the sequence t 15 T2,... ,tm decreases slowly enough, the distribution, 
which starts in the r = 00 Boltzmann distribution, will at step k be close to the Boltzmann 
distribution for temperature rfc. The question is how many transitions must be made in order 
to stay near the Boltzmann distributions as the temperature is lowered down to near 0. If 
the total number of steps required grows only as a polynomial in n, we say that the simulated 
annealing algorithm succeeds. 

4 Symmetrized cost functions 

In this section we discuss cost functions of the form (0) that are functions only of the 
Hamming weight w = Z\ + z 2 + • • • + z n as a result of all bits being treated symmetrically. 
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As an example, let /13 depend on three bits z, z', and z" with the form 



h(z,z',z") 



(0 z + z' + z" = 

q z + z' + z" = 1 

1 z + z' + z" = 2 

1 2 + z' + z" = 3 



(11) 



where q is an integer greater than or equal to 3. Now let 



hfa, ... ,z n ) = ^ h z( 

2!j, Zj, Zfc) 

i<j<k 



(12) 



which gives 



h(w) = h(zi, ... ,z n ) 



\w(n — w)(n — w — 1) + \w{w — 1) (n — w) + — l)(w — 2) 



We can also write 



where 



9W 



h{w) = (^fg{w/n) + 0{n 2 ) 



Aqu{l - uf + 4u 2 (l - u) + f w 3 . 



(13) 



(14) 



(15) 



For n large the form of g(u) determines the behavior of the classical and quantum algorithms. 
In Figure [I] we plot g(u) for q = 3. Note that the global minimum of g is at u = 
corresponding to the unique input Zx, . . . , z n = 0, . . . , 0. The derivative g'{l/2) is negative 
and there is a local minimum at u = 1 corresponding to Z\, . . . , z n = 1, . . . ,1. 

We now look at the quantum algorithm applied to this problem. The cost function ( |i~3"D 
corresponds to the Hamiltonian 



2 V2 



7? 



l(^-s z \[-- + s z \[-- + s. 



S z 



n 



+ -(--s^ 

6 V2 



n 

2 ~ 



S z -1 



where 



n 

2"" 



J \z\, ■ ■ ■ ,z n ) . 



S z \ z li ■ ■ ■ 1 Zn) 

We choose for the initial Hamiltonian (following [|TJ, formula 2.28) 

H(0) 



(16) 



(17) 



77 — 1 

2 



--3 a 
2 
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Figure 1: The leading term in the scaled cost as a function of the scaled Hamming weight u = w/n. 



where 



(0 

X 



(19) 



i=l 



Choosing a linear interpolation between H(0) and Hp gives 

H(s) = (l-s)H(0) + sH P for0<s<l. 



(20) 



Since H(s) depends only on the two operators S x and S z , we start our analysis by considering 
the states |#) defined by 



71 

(sm9S x + cos 9 S z )\9) = -\9) . 



Note that 



77/ 

(e\s x \e) = -sine 



(e\s 2 x \e) = (^mne) 2 + o(n) 

(9\S Z \9) = -cos9 

(9\S 2 Z \9) = (|cos0) 2 + 0(n) 

(9\S 3 z \9) = (^cos9y + 0(n>) 



(21) 

(22a) 
(22b) 
(22c) 
(22d) 
(22e) 



Now define 



V{9, s) = 2(1 - s)(l - sin 9) + sg{\{! - cos 9)) 



(23) 
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where g is given in ([To]) . Using ([I 

V{9, a 



(f20|), and fl22|) we have 



-) 3 (e\H(s)\e) + o(i/n) 



(24) 



where we have also used the operator relation ty = | — S^. The function V(9 } s), for each s, 
represents a large- n "effective potential". For n large, for each s, the ground state of H(s) 
is well approximated by \9 m ) where 9 m minimizes V(9,s). 



s=0 



s=.2 



s=.4 




1 2 3 

theta 

Figure 2: The effective potential versus 6 changes shape as s increases. 

In Figure | we plot V{9, s) versus 9 for s = 0, 0.2, 0.4, 0.4341, 0.46, and 1. At s = the 
minimum is at 9 = n/2 corresponding to Hamming weight n/2. Because ^'(1/2) < 0, this 
minimum moves to 9 bigger than 7r/2 for small s as can be seen in the s = 0.2 panel. As s 
increases, a local minimum forms at values of 9 less than n/2. There is a critical value of s, 
s* = 0.4341, at which the two minima occurring at 9\ < n/2 and at 9 2 > n/2 are degenerate; 
that is, V(9i, s*) = V(92, s*). For s > s* the global minimum is always at a value of 9 less 
than n/2. For s = 1 this global minimum is at 9i = 0. 

For the quantum adiabatic evolution algorithm to work, the quantum state of the system 
should remain in (or very near) the instantaneous ground state of H(s) with s = t/T. 
For s < s* the ground state is smoothly varying. However, at s = s* the ground state 
changes from ~\9 2 ) to For the quantum system to go from |# 2 ) to requires 

quantum tunneling through the barrier separating the minima at 9\ and 9 2 . This takes a 
time exponential in n and implies that the gap of H(s*) is exponentially small. Using a 
combination of standard large spin and instanton methods 0] we can estimate the dominant 
exponentially small term. We obtain 



gap ~ poly(n)e 



-Son 



(25) 
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Figure 3: The log of the gap plotted versus n. The straight line behavior indicates that the gap is 
decreasing exponentially in n. 



where 



So 



p(9)sm9d9 



(26a) 



and 



<p(0) = cosh 



1 + 



v(e, a *)-v{e U 8*) 



2(1 



sin# 



(26b) 



For q = 3, we find s* = 0.4341 and Sq = 0.2015. We can also determine the minimum gap 
of H(s) by numerical diagonalization. Because of the symmetry of H and the initial state, 
although the full Hilbert space is 2 n dimensional, the evolution takes place in the (n + 1)- 
dimensional subspace of total spin n/2. In Figure |^ we show the log of the gap versus n for 
n between 60 and 100. At n — 100 the minimum gap occurs at s — 0.4338. Fitting the gap 
to the functional form An p e~ cn and finding the best values of A, p, and c gives c = 0.2021, 
which certainly supports our approach. For q = 5 we have So = 0.4946 and c = 0.4966, and 
for q = 7 we have S = 0.7008 and c = 0.6930. 

Formula ( f25|) is valid for any example where the cost function depends only on the 
Hamming weight w and we can write 



h(w) = n r g{w/n) + 0{n 



r-l\ 



(27) 



where the global minimum of g(u) occurs at u < 1/2 and g'(l/2) < 0. This guarantees that 
there is a value s* such that V(8, s*) given by (|23|) has degenerate minima at 6\ < ir/2 and 
9 2 > tt/2. Then (EBT) can be used to calculate the coefficient of n in the exponent of 
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We now study the performance of simulated annealing for this kind of problem. Again 
we start with a cost function h(z±, . . . ,z n ), which is a sum of local terms but because of 
symmetry only depends on the Hamming weight, so we can write the cost function as h(w). 
Now for any temperature r, the Boltzmann distribution (|8]) depends only on the Hamming 
weight of the string w = Z\ + z 2 + • • • + z n . The number of strings with Hamming weight w 
is ( ) ■ We can define the (unnormalized) probability 

PrH = (l)ex P {-h(w)/r) (28) 

which gives the probability of finding a string with Hamming weight w at temperature r in 
the Boltzmann distribution. 

We also define the "Free energy", F(w,r), by 

ex V (-F(w,r)/r) = Q exp(-h(w)/r) . (29) 

Note that for any r, the most likely value of the Hamming weight occurs at the value of w 
that minimizes F(w,t). Now we further specialize to the case ([13]) so that (|14]) and ( |T5| ) 
apply. If we rescale the temperature and let a = 8r/n 2 then we can write 

F(w, an 2 /8) = [-) F{u, a) + 0{n 2 ) (30) 

where 

J-(u, a) = aUl — u) log(l — u) + u logw] + g{u) (31) 

and again u = w/n. 

For n large, for any value of the temperature r = an 2 /8, the Boltzmann distribution 
is peaked on those strings with Hamming weight nu where u minimizes (0). At infinite 
temperature, the first term in fl3~l| ) dominates and is minimized at u = 1/2. As a is decreased 
the minimum moves to a value of u above 1/2. In Figure |] we show JF(-u, a) /a versus u for 
a = 100, 10, 3.4, 3.0574, 2.8, and 0.01. The similarity with Figure Q is apparent. The critical 
value of a, a* = 3.0574, is where the global minimum switches from a value of u above 1/2 
to a value below. To follow the Boltzmann distribution, in the simulated annealing Markov 
chain, the Hamming weight, nu, must change from the value corresponding to the larger 
local minimum to the smaller while r is approximately a*n 2 /8. 

This order-n change in the Hamming weight requires exponentially many steps: The sign 
of the derivative -^T{u, a*) determines whether the Markov chain is more likely to move left 
or right. For large n, if in the steps of the annealing process the bits are chosen uniformly 
at random, w is more likely to move to w + 1 than to w — 1 if t^J-{u, a*) < 0. This means 
that climbing the hill separating the local minima is exponentially unlikely in polynomially 
many steps. 

As long as g(u) has a global minimum at u < 1/2 and g'(l/2) < simulated annealing, 
like quantum adiabatic evolution, will require a time exponential in n to succeed. The method 
used in this section, tracking the local minima of the effective potential and the free energy, 
can be used to show the similarity of the performance of quantum adiabatic evolution and 
simulated annealing for many examples with a cost that depends on the Hamming weight 
alone. 
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Figure 4: The scaled free energy J- /a plotted as a function of the scaled Hamming weight changes 
shape as the temperature decreases. Compare to Figure 0. 



5 The Hamming weight with a spike 

Here we consider the cost function 



h(w) 



UK 



n. 



w 7^ n/A 
n/A 



w 



(32) 



where again w is the Hamming weight of n bits and also n is taken to be a multiple of 4. 
Note that this cost function cannot be written as a sum of terms each of which involves 
only a few bits. The global minimum of h(w) is at w = and there is a local minimum at 
w = (n/A) + 1. 

It is easy to see that simulated annealing fails to find the global minimum in polynomial 
time. For r > l/log3, the Boltzmann distribution concentrates at Hamming weights larger 
than n/A. For r < l/log3 the Boltzmann distribution concentrates at Hamming weights 
smaller than n/A. At temperatures ~ 1/ log 3, according to the acceptance rule fllPp, a string 
with Hamming weight (n/A) + 1 will only flip to a string with Hamming weight n/A with a 
probability that is exponentially small in n. Thus a simulated annealing algorithm running 
for only polynomial time gets hung up in the false minimum of h. 

In quantum mechanics, it is well known that a low-energy particle can penetrate a high 
potential barrier. If the barrier is thin, the penetration probability need not be small. This 
turns out to be the case here. As in the examples of Section |], since the cost is a function 
only of the Hamming weight, the quantum evolution takes place in an (n + l)-dimensional 
subspace. The gap can be calculated by explicit diagonalization (using a technique similar 
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Figure 5: The Hamming weight with a spike. The straight line behavior on the log-log plot indicates 
that the gap decreases polynomially in n. 

to the one in U, section 4.2). The result is that 

gap(n) « 1.35 n- 1/2 (33) 

and, for n large, occurs at s = 0.366. We can also find the gap by numerical diagonalization; 
see Figure [5], which plots loggap(n) versus log(n) for 500 < n < 848. The straight line fit 
has slope —0.474, consistent with the predicted —0.5 in (|33|). For n = 848, the gap occurs 
at 0.3675, close to the predicted 0.366. 

6 The bush of implications 

In the previous section we discussed an example with a cost function that treats all bits 
symmetrically and depends only on the Hamming weight of the n bits. Here we describe an 
example where the cost function depends on the Hamming weight of n bits, but also on the 
value of a single additional bit. As in Section we consider the effective potential and free 
energy. In the example in this section, quantum adiabatic evolution succeeds in polynomial 
time while simulated annealing does not. 

We label the n + 1 bits 0, 1, 2, . . . , n. The cost function is given by 

n 

h(z , zi, . . . , z n ) = ^2 ^o(l - Zi) + (1 - z Q ) . (34) 

i=l 

The terms zo(l — Zj) can be viewed as "imply" clauses, since they are unless z = 1 
and Zi = 0. The unique string with h = has zq = 1 and all Zi — 1. If zq — 0, then 
h — 1 no matter what the values of Zi, . . . ,z n are. If Zq — 1, then h equals n — w where 
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1 2 3 ■ ■ • n 
Figure 6: The bush of implications. 



w = Z\ + z 2 + 
as 



and as before 



■ ■ + z n is the Hamming weight of the string zi,... ,z n . We can rewrite (j34|) 



with 



and 



h(z , w) = z (n — w) + l — Zo 



h(z ,w) = ng(z , — ) + 0(1) 



g(0,u) = 



(35) 



(36) 



(37a) 



0(1, tt) = 1 - u . 



(37b) 



Note the degeneracy of the minima of g; this is removed by including the lower order term 
(l-*o). 

As before, to determine the performance of annealing, we calculate the free energy J 7 , 
defined by 



and 



We get 



and 



exp(-F Z0 (w,t)/r) = exp(-h(z ,w)/r) 



F Z0 (w,r)=n^ Z0 i K -,r)+O(l) . 



!Fq{u, t) = r[(l - u) log(l - u) + wlogw] 



Fiiu, t) — t[(1 — u) log(l - u) + ulogu] + (1 — 



u 



(38) 



(39) 



(40a) 



(40b) 



In Figure we plot T Zo {u,T)jr, r = 100,5, 1, and 0.05. We see that with the 0(1) terms 
neglected, the value z = 0, w = | minimizes the free energy for all r, < r < oo. 
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Figure 7: The bush of implications. The scaled free energy plotted versus the scaled Hamming 
weight as the temperature changes. In each panel the top curve corresponds to zq = 1 and the 
lower corresponds to zq = 0. The global minimum never moves from zq = 0, u = 1/2. 

Can the addition of the correction, which adds 1 to Fq(w,t) for all (w,t), make any 
difference in a simulated annealing process? For r = 0, the 1 — zq term means that the 
global minimum of h occurs at Zq = 1, w = n. For r ~ 1, a simulated annealing process 
that is following the Boltzmann distribution will be concentrated at Zq = 0, w = | . There 
are two routes from zq = 0, w = ~ to z = 1, w = n. One route is to jump from z = to 
Zq = 1 while w ~ 77. But the probability of accepting this change is exponentially small in 
n for t < 1, so this route is exponentially unlikely in polynomially many steps. (Note that 
changing the rate at which bit is nipped does not help.) The other route is to maintain 
zq = while w increases. But this is exactly the kind of hill climbing in free energy [this time 
JF (u, r)] that cannot occur in polynomial time as in the example in the previous section. In 
conclusion, simulated annealing fails to reach the true minimum of h at Zq — 1, w — n in 
polynomial time. 

We now turn to the behavior of a quantum adiabatic algorithm for the bush of implica- 
tions. Using spin operators we write the quantum Hamiltonian corresponding to h of ( ^4[) 

as 

Hp = \{l-^){l + S z ) + \{l + of) (41) 

where S z is the z component of the total spin for bits 1 through n, whereas ^<j[°^ is the spin 
component of the spin of bit 0. Now bits 1 through n are each involved in one term in ( |34"|) 
whereas bit is in n + 1 terms. Accordingly, we get for H(0) [see |L|, formula (2.22)] 

H(0) = ( n +l)l(i -,(«>) + (| - Sa ) . (42) 
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Figure 8: The bush of implications. The global minimum of the effective potential plotted as a 
function of s. 



For the interpolating Hamiltonian we take 

H(s) = (1 - s)H(0) + sH P . (43) 
Consider the states a, 0) = a)<g> 0) where a) is the bit state 

\a) = (sma/2)\z = 0) + (cosa/2)|z = l) (44) 
and |0) is made of bits 1 through n and is defined by (pip . For the effective potential we set 

ol 9 

V(a, 0, s ) = (1 - s) [2 - sin a - sin 0] + 2s cos 2 - cos 2 - (45) 

which gives 

V(a,9,s) = -(a,9\H(s)\a,9) + o(— ) . (46) 

Note that V(a, 0, s) = V(9, a, s) for all s. 

In Figure || we plot (a(s),9(s)), the global minimum of V(a,9,s), as a parameterized 
curve for < s < 1. At s = the global minimum occurs at a = n/2, 9 = tt/2, which 
corresponds to a uniform superposition for bit and a Hamming weight concentrated at n/2 
for bits 1 through n. As s increases the minimum is at a(s) = 9(s) until a critical value 
s = s* = 0.805. For s > s* there are two global minima symmetric under interchange of a 
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Figure 9: The bush of implications. The difference between the two lowest eigenvalues of the 
Hamiltonian as a function of s, when n = 300. Note that the minimum occurs near the predicted 
point s* = 0.805. 



and 9. We go beyond the effective potential and analyze the quantum corrections for s = s*, 
using a representation in terms of a 2-component (to represent the state of bit 0) anharmonic 
oscillator (to represent the fluctuations about the state We find an expansion in 

powers of n" 1 / 3 . The symmetry under interchange of a and 9 is broken by the quantum 
corrections, and the branch heading towards a = rr/2, 9 = tt has the lower energy. The 
adiabatic evolution is therefore heading towards the state that is a uniform superposition for 
bit and has Hamming weight n for bits 1 through n. This analysis so far has neglected the 
|(1 + o"i°' ) ) term in Hp because it is lower order in n. For (1 — s) of order 1/n, the effect of 

this term in H(s) is of the same order as H(0). The net effect of this term is to rotate bit 
to |^o — l) m an n-independent fashion. 

According to this analysis the minimum gap appears at s = s* + 0(n~ 2 / 3 ), where s* is 
the bifurcation point in Figure [8l We calculate 



gap(n) 



0.5782 



B 



n 



n 



2/3 



+ o(n- 1 ; 



(47) 



We can confirm this analysis numerically because the Hamiltonian H(s) acts in a reduced 
Hilbert space of dimension 2{n + 1). In Figure ^| we show the difference between the lowest 
two eigenvalues E\(s) — Eq(s) for 0<s<latn = 300. Note that the minimum difference 
occurs at s = 0.800. In Figure [10] we plot n 1 ^ 3 ■ gap(n) versus n _1//3 for 300 < n < 1540. 
According to prediction (f47D this should be a straight line (with corrections of order n~ 2 / 3 ) 
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Figure 10: The bush of implications. The gap multiplied by n 1//3 and plotted versus n~ 1//3 . The 
intercept of a straight line fit is 0.5812, compared to the predicted 0.5782. 



with an intercept of 0.5782. A straight-line fit to the data in Figure [10] has an intercept of 
0.5812. 

It is worth noting that the choice of H(0) in (f42|) was made so that H(s) would be the 
sum of terms 



Ms) 



( 4 * 



c 



where each H c (s) depends only on clause C (see [I], formula 2.28). Suppose instead, for the 
bush of implications, we write 



H(0) = X(n + 1)1(1 + 



(49) 



with A a parameter. We can recalculate the effective potential and draw the picture analogous 
to Figure §. If A > 1, the gap for n large is a positive constant. If A < 1, the gap is 
exponentially small. Note that the choice A = l/(ra+l) corresponds to a magnetic field in 
the ^-direction with the same strength for each bit. This case, mentioned in 0, is an example 
of a local problem with a choice for a quantum adiabatic algorithm that requires exponential 
time. This illustrates that the cost function alone does not determine the performance of 
the quantum algorithm but rather that the choice of the path to Hp is sometimes crucial. 
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7 Conclusion 

In this paper we first compared the behavior of quantum adiabatic evolution and classical 
simulated annealing in searching for the global minima of cost functions that depend on 
n bits symmetrically. When the cost depends smoothly on the Hamming weight of the n 
bits, quantum adiabatic evolution and simulated annealing typically perform similarly. 

We have also shown, however, that there are examples where this similarity breaks down. 
In the spike example, the cost function has a barrier that can be penetrated quantum me- 
chanically but not classically, so the quantum algorithm succeeds in polynomial time whereas 
annealing does not. In the bush example, the addition of a single spin-1/2 leads to quantum 
behavior with no classical analogue. These examples suffice to show that there is no theorem 
that says quantum adiabatic evolution will fail when simulated annealing fails. It remains an 
open question whether they also indicate that quantum adiabatic evolution algorithms will 
succeed on problems of computational interest that cannot be solved efficiently by classical 
local search. 
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